Nontrivial in-plane-magnetic-field dependence of THz wave emission 
from intrinsic Josephson junctions controlled by surface impedance 
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Recently THz wave emission from intrinsic Josephson junctions (IJJs) was confirmed without 
external magnetic fields, and the surface impedance Z was found out to play a crucial role. When in- 
plane magnetic fields are introduced, field dependence of emission intensity also strongly depends on 
Z. For Z = 1, emission becomes the strongest when the wavelength of electromagnetic wave in IJJs 
coincides with the distance between Josephson vortices (JVs). For Z > cavity resonance modes 
are stabilized and the fundamental mode gives the strongest emission. There occurs a dynamical 
phase transition between the 7r-phase-kink state (stable without external fields) and in-phase state 
as the in-plane magnetic field increases. Although the incommensurate-phase-kink state between 
these two states would be an artifact for small number of layers used in simulations, this artificial 
state almost vanishes for Z ~ 50 ~ 60, where field profile of emission intensity changes from a 
characteristic peak around 3/4 JVs per layer for smaller Z to monotonic decrease for larger Z. This 
crossover induced by Z may explain recent controversial experiments on the field profile of emission. 

PACS numbers: 74.50.-|-r, 85.25.Cp, 74.25.Nf 



Introduction. Although THz wave emission from in- 
trinsic Josephson junctions (IJJs) had been investigated 
in in-plane external magnetic fields, experimental real- 
ization in such a condition [Ij had been quite difficult. 
Recently evident THz wave emission was observed with- 
out external magnetic fields. 0, Q In these experiments 
about a thousand of layers were included in IJJs so that 
the velocity of Josephson plasma mode is already as fast 
as that of light in IJJs, and only the fundamental (n = 1) 
cavity resonance mode has been observed so far. Then, 
possible emission states were investigated theoretically, 
and the uniform in-phase state 0, [Hjl and the 7r-phase- 
kink states with symmetry breaking along the c axis 0] 
have been proposed. The present author showed Q that 
these states are both stationary according to the bias 
current J and surface impedance Z , and drew the dy- 
namical phase diagram in the Z-J plane. Quite recently 
emission in in-plane fields has been investigated again. In 
an experiment intensity decreases monotonically as the 
in-plane field increases, [9] while in another experiment 
the emission intensity seems to have some characteristic 
peaks in the field profile. 10] The present study suggests 
that these two experiments may not be controversial. 

Model and formulation. As long as thermal fluctu- 
ations are not taken into account in the modeling, di- 
mensional reduction along the in-plane magnetic field is 
justified. When the direction of the field is chosen as the 
y axis, the basic equations are given as follows, jTlj 



'1+1,1 



(i-CA(2))(a,,i?,Vi,+ 

+ sin ipi+i,i 
(l-aA(2))£;; 



(1) 
(2) 



where the subscript "I + 1,1" denotes quantities in the 
insulating layer between the Z-th and {I + l)-th super- 
conducting layers, and the operator A*^^^ is defined in 



A(2)x,+i,i EE Xi+2,i+i - 2Xi+i^i + Xij^i. Definitions of 
scaled quantities are given in Eqs. (7)-(10) of Ref. Q- For 
example, length scale is scaled by the penetration depth 
in x' , time scale by the inverse of plasma frequency in t' , 
and bias current by the critical current in J'. In addition 
to the scaled electric field -Ej^^i ; and the gauge-invariant 
phase difference ipi+i,i, the scaled magnetic field is ob- 
tained from dx''4'i+i,i = (1 — C^'^')5i+i Using the ma- 
terial parameters of Bi2Sr2CaCu208 given in Ref. fll|. 
we have a large inductive coupling ^ = 4.4 x 10^ and a 
small capacitive coupling a = 0.1, and /3 ~ 0.02 is taken. 

Width of the junction = 86/im is comparable to 
the size of samples in experiments, and for this size 
By ~ 0.02T corresponds to the condition "one JV per 
layer" . Although a thousand of layer is essential for ob- 
taining the plasma velocity as fast as that of light in 
IJJs, it is almost impossible to arrive at stationary states 
for such a large number of layers numerically at present. 
Then, the periodic boundary condition along the c axis is 
introduced instead. In such a condition plasma velocity 
of the stationary state automatically coincides with that 
of light in IJJs, though vanishing amplitude of electro- 
magnetic waves on surfaces is totally neglected. Even if 
a thousand of layers are stacked, thickness of IJJ is 
much smaller than the wavelength of emitted THz wave 
A. In such a case emission is given by a point-like source 
(after two-dimensional reduction, by a line-like source in 
original three dimensions) and weaker than that from in- 
finite sources. Such effect is approximately included as: 
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with the dielectric constants of IJJs e' 



and dielectrics e^. 
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FIG. 1: Field dependence of maximum intensity in the n — 1 
cavity resonance mode for = 4 (circles) at the optimal bias 
current (squares). Dips of the curves stand for the dynamical 
phase transitions between the (a) 7r-phase-kink (tt-PK), (b) 
incommensurate- phase-kink (IPK), and (c) in-phase states. 
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FIG. 2: Snapshots of gauge-invariant phase differences in all 
layers at (a) By = 0.0033T, J = 0.166Jc (in the tt-PK state 
near the phase boundary; two ±7r kinks), (b) By = 0.006T, 
J — 0.159Jc (in the IPK state near the phase boundary; three 
-|-(l/2)7r kinks and one -(3/2)7r kink), and (c) By = 0.015T, 
J — 0.167Jc (in the in-phase state near the emission peak). 



respectively. Here the dynamical parts of scaled bound- 
ary magnetic field B'l^i ; and electric field E[^^ ^ are 
related with each other by z « A/L^. [12] Although 
this naive evaluation of z might be modified by excess 
magnetic fields from vertical directions omitted in two- 
dimensional modeling, effect of impedance mismatch 
on the edges is also included in Z. Width of the sample 
Lx is divided into 80 numerical grids, and numerical cal- 
culations are based on the RADAU5 ODE solver. 

Numerical results for Z — 1. This condition is typi- 
cally realized for z — 1 and — ^'d- Dependence of emis- 
sion intensity on in-plane magnetic field and bias current 
was already investigated in Ref. [15]. Here the results are 
summarized for comparison with those for larger Z. 

Aside from the recoupling region for small current, 
emission intensity takes the maximum value at the point 
where the wavelength of electromagnetic wave in IJJs 
and the distance between JVs coincide. From this con- 
dition and the ac Josephson relation, the current at the 
maximum intensity is proportional to the external mag- 
netic field. The maximum intensity is independent of the 
current until it reaches the critical one Jc. On the con- 



trary to previous numerical studies for Z = I, ll|, [16[ no 
symptom of the cavity resonance remains in this simu- 
lation. In the previous studies number of layers is finite 
and effect of surface electrodes (penetration of supercon- 
ductivity into electrodes) is taken into account, and this 
effect played a role of "cavity" even for Z = 1. Such an 
effect does not exist under the periodic boundary condi- 
tion along the c axis. The frequency spectrum is quite 
sharp at the peak of emission intensity for each value of 
the in-plane magnetic field, and the frequency at emission 
peak is proportional to the in-plane magnetic field. 
Numerical results for Z = 10. Similarly to the emission 



without external magnetic fields, the emission in the in- 
plane magnetic fields sharply depends on the value of 
the surface impedance Z. The cavity resonance modes 
are stable for Z > 3, and the fundamental {n = 1) mode 
becomes the strongest in in-plane magnetic fields. In the 
present article only the n = 1 mode is considered. 

Field dependence of the maximum intensity with the 
optimal bias current for Z = 10 is shown in Fig. [T] by 
circles. Two dips in this figure divide the three dy- 
namical states characterized by the structure of gauge- 
invariant phase difference shown in Fig. [5] (a) the tt- 
phase-kink (tt-PK) state stable without external fields, 
(b) the Incommensurate-phase-kink (IPK) state in the in- 
termediate field region, and (c) the in-phase state. Field 
dependence of the current deriving the maximum inten- 
sity (optimal current) is also plotted in Fig.[T]by squares. 
Its shape is quite similar to that of the maximum inten- 
sity, and this field profile will not be shown hereafter. 

The TT-PK state is expected to be stable even for larger 
number of layers because stability of this state has al- 
ready been confirmed without magnetic fields. If the in- 
phase state is stable, its property is expected to be inde- 
pendent of the number of layers. Then, peculiar layer- 
number dependence of physical properties may be seen 
in the IPK state. Field dependence of the maximum 
intensity for = 4, 6, 8 and 12 is displayed around 
this region in Fig. [31 and typical snapshots of the gauge- 
invariant phase difference are shown in Figs. |l{a), IH^b) 
and[4{c) for A^ = 6, 8 and 12, respectively. These exam- 
ples suggest that the IPK state for A(= even) layers is 
characterized by ±(1 -I- 2/N)'k- and =f(1 — 2/A^)7r-phase 
kinks with numbers (A^/2 — 1) and (7V/2-I-1), respectively. 

Layer-number dependence of the phase boundaries 
shown in Fig. [3] is not trivial. The boundaries are de- 
termined by two factors: stability of each state and dif- 
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FIG. 3: Field dependence of maximum intensity in the n — 1 
cavity resonance mode for A'^ = 4 (open circles), 6 (full trian- 
gles), 8 (full squares), and 12 (full circles). Dash-dotted and 
dashed lines denote boundaries between the vr-PK and IPK 
states and between the IPK and in-phase states, respectively. 



ference between the states. The tt-PK and in-phase states 
are well-defined independent of the number of layers N, 
while the IPK state varies as does. Fig. [3] indicates 
that intensity of emission monotonically decreases as N 
increases in the IPK state, which suggests that this state 
becomes more and more unstable as N increases. It 
would be why the boundary between the IPK and in- 
phase states shifts to lower magnetic fields for increasing 
N. The boundary between the tt-PK and IPK states is 
more nontrivial because similarity between the tt-PK and 
IPK states also becomes stronger as N increases, and fi- 
nally arrives at the same state in the iV — > cxd limit. 

Then, the boundary fields between the tt-PK, IPK and 
in-phase states are plotted versus inverse of layer num- 
bers in Fig. Hljd). Since the boundary field between the 
TT-PK and IPK states is expected to converge to By — 
in the N ^ oo limit, layer-number dependence of phys- 
ical quantities can be checked by this condition. Fig. 
mjd) shows that the naive plot is not bad (dash- 
dotted line versus full circles corresponding to the data 
for = 6, 8 and 12). When the data for the boundary 
between the IPK and in-phase states (denoted by open 
circles) are extrapolated to iV — > cxd assuming similar 
N dependence (dashed lines), they clearly converge to a 
finite value. In summary, when the small-size effect of 
layer numbers is removed, there occurs one dynamical 
phase transition between the tt-PK and in-phase states, 
and in this tt-PK state the alternate arrangement of the 
-|-TT- and — TT-phase kinks [7] would be observed only in 
the vicinity of zero external field, and there would be a 
crossover to random arrangements of these two kinks. 

Numerical results for other Z . As shown so far, emis- 
sion behaviors ioi Z — 1 and 10 are quite different. This 
situation is similar to the zero-field case, and cavity- 
resonance properties become apparent for Z > 3. Field 
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FIG. 4: Typical snapshots of gauge-invariant phase differences 
in all layers in the IPK region for (a) A'' = 6 (four +(2/3)Tr 
kinks and two -(4/3)Tr kinks), (b) iV = 8 (five +(3/4)Tr kinks 
and three — (5/4)Tr kinks), and (c) TV = 12 (seven -|-(5/6)Tr 
kinks and five — (7/6)Tr kinks), (d) Boundary fields between 
the TT-PK and IPK states (full circles) and the IPK and in- 
phase states (open circles) versus inverse of layer number A''. 
By naive linear extrapolation to the N ^ oo limit, the for- 
mer (dash-dotted lines) and the latter (dashed lines) seem to 
converge to -By = and a finite value of By, respectively. 



dependence of maximum intensity of emission in this re- 
gion of the n — 1 mode for A^ = 4 is given in Fig. [DJa) 
for Z = 3, 10, 30, 100 and 300 together with transition 
fields for each Z. This figure clearly indicates large Z 
dependence of the field profile of maximum intensity of 
emission. Although emission peak in the IPK region for 
Z — 10 and 30 is an artifact of small number of layers 
(A^ = 4), there is another emission peak in the in-phase 
state at By w 0.015T for Z = 10 and 30. On the other 
hand, there are no characteristic peaks for Z = 100 and 
300; By = 0.02T (corresponding to 1 JV per layer) still 
belongs to the IPK region and even the field profile be- 
comes almost field independent for Z = 300. 
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FIG. 5: Field dependence of maximum intensity for A'^ = 4 
for various values of the surface impedance Z: (a) for Z = 3 
(stars), 10 (open circles), 30 (full circles), 100 (full squares) 
and 300 (full triangles), (b) expanded figure for Z — 40 (open 
triangles), 50 (full squares), 60 (open squares) and 70 (full 
triangles) . In each figure dash-dotted and dashed lines denote 
boundaries between the vr-PK and IPK states and between the 
IPK and in-phase states, respectively. 



consequence of small difference in Z in the vicinity of the 
crossover point. This argument also suggests that val- 
ues of Z in experiments may be much smaller than that 
naively expected from the relation z « X/Lz ov Z Ki 500. 
It means that the argument by Tachiki et al. [l3| that Z 
is effectively reduced by magnetic fields from vertical di- 
rections may be partially true, but that such effect is not 
so strong as they considered; they argued that Z may be 
reduced so much that the in-phase emission is observed. 

Summary. THz wave emission from intrinsic Joseph- 
son junctions in in-plane magnetic fields is investigated 
numerically. Frequency-tunable emission by in-plane 
fields is observed for small surface impedance Z < 3, 
and cavity-resonant emission is observed for Z > 3. 
Experimentally-accessible n ~ 1 cavity resonance mode 
is the strongest, and there occur field-induced dynami- 
cal phase transitions between the 7r-phase-kink (tt-PK), 
incommensurate-phase-kink (IPK) and in-phase states as 
the in-plane magnetic field increases. Investigation on 
size dependence of physical quantities on layer numbers 
N shows that IPK state may be an artifact of small N 
{N w 1000 in experiments) and that only the tt-PK and 
in-phase states may be stable for large N. In the field 
profile of maximum intensity there exists a characteristic 
peak around 3/4 JVs per layer [l3| for Z < 40, while it 
decreases monotonically [oj for Z > 70. This crossover 
of the field profile is characterized by sharp shrink of the 
IPK region even for iV = 4 for Z « 50 ^ 60, which 
suggests that controversy of the field profile in recent ex- 
periments [loj may be explained by slight difference in 
the surface impedance due to experimental conditions. 
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Then, Z dependence of the field profile of maximum 
intensity between Z = 30 and 100 is analyzed precisely. 
The field profile for Z = 40, 50, 60 and 70 is displayed 
in Fig. [S{b) together with transition fields for each Z. 
While the characteristic peak at By w 0.015T is still 
visible for Z = 40 as for Z = 30, it almost disappears for 
Z = 50 and 60, where the IPK region shrinks only around 
By = 0.012T. For Z = 70 the IPK region broadens and 
even By = 0.02T is included in this region as for Z = 100. 
That is, there exists a crossover between two different 
field profiles of maximum intensity at Z « 50 ~ 60. 

Discussions. Quite recently there appear two experi- 
ments for THz wave emission from IJJs in in-plane fields 
based on the same samples for emission without mag- 
netic fields. Welp et al. §| observed monotonic decrease 
of emission intensity as the field increases, while Yamaki 
et al. pii found characteristic peaks in the field profile of 
emission intensity. If such experimental conditions cor- 
respond to the crossover region at Z « 50 ^ 60, these 
two experiments may not be a serious controversy but a 
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